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Abstract 

Using data collected on the ip(3770) resonance and near the D*^D^ peak production en- 
ergy by the CLEO-c detector, we study the decays of the possible D — ► PP modes and report 
measurements of or upper limits on all branching fractions for Cabibbo-favored, singly-Cabibbo- 
suppressed, and doubly-Cabibbo-suppressed D — > PP decays except modes involving (and 
except D° — > K + ir~). We normalize with respect to the Cabibbo-favored D modes, D° — > K~ir + , 
D+ -► K-7T+7T+, and -► 
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I. INTRODUCTION 

There are many possible exclusive decays of charmed D mesons to a pair of mesons from 
the lowest-lying pseudoscalar meson nonet. The decay can be to any pair of K + , K~, n + , n~, 
T), r]', 7T°, K°, or K°, with total charge or ±1. Measurements of the complete set of decays 
can be used to test flavor topology and SU(3) predictions and to specify strong phases of 
decay amplitudes through triangle relations [1]. Moreover, many CP asymmetries (expected 
to be less than O(l0~ 3 ) in the Standard Model) can be studied. The detectable neutral 
kaons are K° s and K\, not K° and K°, so the observable decays are XK S and XK\. In this 
study, we consider only Kg, not K\, and report all branching fractions for Cabibbo-favored, 
singly-Cabibbo-suppressed, and doubly-Cabibbo-suppressed D — > PP decays except modes 
involving K\ and except the doubly-Cabibbo-suppressed decay D° — > K + n~ . We normalize 
with respect to the Cabibbo-favored D modes, D° — > K~n + [2], D + — > K~ir + n + [2], and 
Df — > K + K s [3]. (More precisely, we normalize the D° — > PP decays with respect to the 
sum of the Cabibbo-favored mode D° — > K~n + and the doubly-Cabibbo-suppressed mode 
D° — > K + 7i~ . The latter is 0.4% of the former.) 

II. THE DETECTOR 

Data for this analysis were taken at the Cornell Electron Storage Ring (CESR) using 
the CLEO-c general-purpose solenoidal detector, which is described in detail elsewhere [4- 
7]. The charged particle tracking system covers a solid angle of 93% of An and consists of a 
small-radius, six-layer, low-mass, stereo wire drift chamber, concentric with, and surrounded 
by, a 47-layer cylindrical central drift chamber. The chambers operate in a 1.0 T magnetic 
field. The root-mean-square (rms) momentum resolution achieved with the tracking sys- 
tem is approximately 0.6% at p — 1 GeV/c for tracks that traverse all layers of the drift 
chamber. Photons are detected in an electromagnetic calorimeter consisting of 7800 cesium 
iodide crystals and covering 95% of 4tt, which achieves a photon energy resolution of 2.2% 
at E 1 — 1 GeV and 6% at 100 MeV. We utilize two particle identification (PID) devices to 
separate charged kaons from pions: the central drift chamber, which provides measurements 
of ionization energy loss (dE/dx), and, surrounding this drift chamber, a cylindrical ring- 
imaging Cherenkov (RICH) detector, whose active solid angle is 80% of 4n. The combined 
PID system has a pion or kaon efficiency > 85% and a probability of pions faking kaons 
(or vice versa) < 5% [2]. The response of the CLEO-c detector is studied with a detailed 
GEANT-based [? ] Monte Carlo (MC) simulation, with initial particle trajectories gener- 
ated by EvtGen [8] and final state radiation produced by PHOTOS [9]. Simulated events 
are reconstructed and selected for analysis with the reconstruction programs and selection 
criteria used for data. 

III. THE DATA SAMPLE 

For D° and D + meson decays, we utilize a total integrated luminosity of 818 pb" 1 of 
e + e~ data collected at center-of-mass (CM) energies near E CTa = 3774 MeV. The data 
sample contains about 2.4 x 10 6 D + D~ events (events of interest), three million D°D° 
events (events of interest), fifteen million e + e~ — > uu, dd, or ss continuum events, three 
million e + e~ — > r + r~ events, and three million e + e~ — > •yip' radiative return events (sources 
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of background), as well as Bhabha events, /z-pair events, and 77 events (useful for luminosity 
determination and resolution studies). For the Df meson decays, we use a data sample of 
e + e- -> D*±DJ events collected at the CM energy 4170 MeV, near D^DJ peak production 
of ~1 nb [10]. The data sample consists of an integrated luminosity of 586 pb" 1 containing 
about 5.4 x 10 5 D^DJ pairs. Other charm production totals ~7 nb [10], and the underlying 
light-quark "continuum" is about 12 nb. Through this paper, charge conjugate modes are 
implicitly assumed, unless otherwise noted. 

IV. PROCEDURE 

A. D° and D+ 

Here we employ a single-tag (ST) technique extensively used by CLEO-c [2, 3, 11, 12], 
pioneered by the Mark III Collaboration at SPEAR for measuring D° and D + branching 
fractions [13, 14], which exploits a feature of near-threshold production of charmed mesons, 
i.e. M bc and AE, see below. 

We formed D and D candidates in all D — > PP decay modes from combinations of n ± , 
K ± , 7T°, Kg, 77, and 7/ candidates selected using the standardized requirements which are 
common to many CLEO-c analyses involving D decays. The ?/>(3770) resonance is below the 
kinematic threshold for DDir production, so the events of interest, e + e~ — > ^(3770) — > DD, 
have D mesons with energy equal to the beam energy. Two variables reflecting energy 
and momentum conservation are used to identify valid D candidates. They are AE = 
T,iEi ~ E heam , and M bc = \] El c&m - (^ Pi) 2 , where E h are the energy and momentum 
of the decay products of a D candidate. For a correct combination of particles, AE will 
be consistent with zero, and the beam-constrained mass M^ c will be consistent with the D 
mass. Candidates are rejected if they fail mode-dependent AE requirements. If there is 
more than one candidate in a particular D or D decay mode, we choose the candidate with 
the smallest |Ai?|. 

B. D+ 

Unlike DD threshold events, conventional AE and M hc variables are no longer good 
variables for D s from D* + D~ decays, as the D s can either be a primary or secondary 
(from a D* decay), with different momentum. We use the reconstructed invariant mass of 
the D s candidate, M(D S ), and the mass recoiling against the D s candidate, M reco [\(D s ) = 

y/ {Eq — Ed s ) 2 — (p — PdJ 2 ) as our primary kinematic variables to select a D s candidate. 
Here (E , p ) is the net four-momentum of the e + e~ system, taking the finite beam crossing 
angle into account, p/j s is the momentum of the D s candidate, E Ds = ^ m 2 Ds + pj- )g , and 
m Ds is the known D s mass [15]. We make no requirements on the decay of the other D s in 
the event. 

There are two components in the recoil mass distribution, a peak around the D* mass 
if the candidate is due to the primary D s and a rectangular shaped distribution if the 
candidate is due to the secondary D s from a D* decay. The edges of M recoi i(D s ) from the 
secondary D s are kinematically determined (as a function of y/s and known masses), and 
aty/s = 4170 MeV, AM Iccoil (D s ) = M Iccoil (D s ) - m D * is in the range [-54, 57] MeV. Initial 
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state radiation causes a tail on the high side, above 57 MeV. We select D s candidates within 
the —55 MeV < AM rccoi i(D s ) < +55 MeV range. This window allows both primary and 
secondary D s candidates to be selected. 

We also require a photon consistent with coming from D* s + — > D+j decay, by 
looking at the mass recoiling against the D s candidate plus 7 system, M recoi i(D s 7) = 

yj (E — E Da — E y ) 2 — (p — p Ds — p 7 ) 2 . For correct combinations, this recoil mass peaks 
at rriD s , regardless of whether the candidate is due to a primary or a secondary D s . We 
require \M recoi \{D s ^) — m Da \ < 30 MeV. This requirement improves the signal to noise ra- 
tio, important for the suppressed modes. Every event is allowed to contribute a maximum 
of one D s candidate per mode and charge. If there are multiple candidates, the one with 
M reco n(D s ^) closest to m Da is chosen. 

C. Common 

Our standard final-state particle selection requirements are described in detail else- 
where [2]. Charged tracks produced in the D decay are required to satisfy criteria based on 
the track fit quality, and angles 6 with respect to the beam line, satisfying | cos6>| < 0.93. 
Momenta of charged particles utilized in D° and D + candidate reconstructions must be 
above 50 MeV/c, while those for D s must be above 100 MeV/c to eliminate the soft pions 
from D*D* and D*D decays (through D* — > ttD). Tracks must also be consistent with 
their coming from the interaction point in three dimensions. Pion and kaon candidates are 
required to have dE/dx measurements within three standard deviations (3a) of the expected 
value. For tracks with momenta greater than 700 MeV/c, RICH information, if available, is 
combined with dE/dx. 

The K s candidates are selected from pairs of oppositely-charged and vertex-constrained 
tracks having invariant mass within 7.5 MeV, or roughly 3a, of the known K s mass [15]. We 
identify n° candidates via n° — > 77, detecting the photons in the Csl calorimeter. To avoid 
having both photons in a region of poorer energy resolution, we require that at least one of 
the photons be in the "good barrel" region, | cos# 7 | < 0.80. We require that a calorimeter 
cluster has a measured energy above 30 MeV, has a lateral distribution consistent with that 
from photons, and not be matched to any charged track. The invariant mass of the photon 
pair is required to be within 3a (a ~ 6 MeV) of the known tt° mass. A ir° mass constraint is 
imposed when n° candidates are used in further reconstruction. We reconstruct rj candidates 
in the decay of rj — > 77. Candidates are formed using a similar procedure as for n° except 
that a ~ 12 MeV. We reconstruct rj' candidates in the decay mode 7/ — > -n +r n^r]. We require 
Im^+Tr-^ — m v /\ < 10 MeV. 

V. RESULTS 
A. D° and D+ 

The Mbc distributions for the D° and D + candidate combinations are shown in Figs. 1 
and 2, respectively. The points show the data and the lines are fits. The normalization 
modes D° — > K~n + and D + — > K~-n + it + are essentially background-free. The backgrounds 
of all modes are well described by the distributions obtained from the AE sidebands. We 
perform a binned maximum likelihood fit to extract the D° or D + signal yield from each 
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FIG. 1: Mb c distributions of D° modes. For each distribution, the points are obtained from the 
AE signal region, the shaded histogram is from the AE sidebands, and the line is the fit. 



Mb c distribution. For the signal, we use an inverted Crystal Ball line shape [16], which 
is a Gaussian with a high-side tail. For the background, we use an ARGUS function [17], 
with the shape parameter determined from the AE sideband Mb c distribution, the high- 
end cutoff given by -Ebeam, and the normalization determined from the fit to the AE signal 
region. Results of the fits are shown in Table I. Table I also includes the detection efficiency 
for each mode. The efficiencies include sub-mode branching fractions [15] and have been 
corrected to include four known small differences between data and Monte Carlo simulation, 
in particular 7r°-finding efficiency 0.96, 77-finding efficiency 0.935, ^ particle identification 
0.995, and particle identification 0.99, data efficiency being smaller than MC efficiency 
by those ratios. 



B. D+ 

The resulting M(D S ) distributions for D s modes are shown in Fig. 3. The points show 
the data and the lines are fits. We perform binned maximum likelihood fits to extract signal 
yields from the M(D S ) distributions. For the signal, we use the sum of two Gaussians for 
the line shape. For the background, we use a second-degree polynomial function. Results of 
the fits and detection efficiencies are given in Table I. 
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FIG. 3: M{D S ) distributions for D s modes. For each distribution, the points are the data and 
the superimposed line is the fit (the dotted line is the fitted background). The distribution for 
— > 7r + 7r° has tighter requirements than the other modes - see text. 
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TABLE I: Observed yields from data and reconstruction efficiencies and their statistical uncertain- 
ties. The efficiencies include sub-mode branching fractions [15] and have been corrected to include 
several known small differences between data and Monte Carlo simulation. 



Mode 




Efficiency (%) 


Yield 


D 




K + K 


57.35 ± 0.16 


13782 ± 136 


U 




K S K S 


22.73 ± 0.13 


215 ± 23 


no 
JJ 






72.68 ± 0.14 


6210 ± 93 


7-1O 




7T 7T 


32. 95 ± U.14 


156/ ± 54 


7-1O 
D 




K 7T + 


65.11 ± 0.15 


150259 ± 420 


nO 

D 




TsO 

K s* 


28.57 ± 0.14 


20045 ± 165 


7-1O 
D 






10.08 ± 0.05 


2864 ± 65 


7-1O 




n 

Tr u r] 


11.97 ± 0.05 


481 ± 40 


JJ 




K%rf 


2.35 ± 0.02 


1 - tit 1 1 ^ r\ 

1321 ± 42 


7-1O 




/ 

7T 77 


2.97 ± 0.02 


159 ± 19 


JJ 




r/rj 


4.35 ± 0.02 


430 ± 29 


u 




1 

r\r\ 


1.06 ± 0.01 


66 ± 15 








54.92 ± 0.16 


231058 ± 515 




~^ 


K S K+ 


36.62 ± 0.15 


5161 ± 86 






+ 


48.69 ± 0.15 


2649 ± 76 


JJ^ 


* 




42.54 ± 0.16 


30095 ± 191 


JJ^ 


* 




A O I nir 

43.29 ± 0.15 


O A O 1 O T 

343 ± 37 








15.95 ± 0.06 


/^r\ 1 A 

60 ± 24 


JJ^ 




4 


18.07 ± 0.06 


2940 ± 68 


n+ 
JJ 




TS+ 1 


4.29 ± 0.02 


23 ± 18 


D+ 




7T 77 


4.ol ± U.U2 


103/ ± 35 


Df 






24.73 ± 0.14 


4076 ± 71 


Df 




7T+7T 


16.60 ± 0.12 


19 ± 28 


Df 






28.15 ± 0.14 


393 ± 33 


Df 






29.57 ± 0.14 


202 ± 70 


Df 






11.40 ± 0.05 


222 ± 41 


Df 




■K + rj 


12.70 ± 0.06 


2587 ± 89 


Df 




K+rj' 


2.87 ± 0.02 


56 ± 17 


Df 




n + rf 


3.28 ± 0.02 


1436 ± 47 



C. Upper Limits 

For most of the D — > PP modes, very clear signals are found in data. We find no 
significant evidence for D + — > K + rj, D + — > K + r/', and Df — > 7r + 7r° decays, and therefore 
set upper limits on their branching fractions. The M bc distributions of D + — > and 
D + — > modes are shown in Fig. 2. Monte Carlo studies indicate that tightening the 

requirements on M reco ii(D s ) to ±10 MeV and M reco ii(D s 7) to ±20 MeV should improve the 
upper limit on Df — > 7r + 7r° decay. Consequently, for Df — > 7r + 7r° (and only — > 7r + 7r ), 
we have applied these tighter requirements. The invariant mass distribution for Df — > 7r + 7r° 
shown in Fig. 3 and the efficiency given in Table I have these tighter requirements. 
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D. Background from Non-resonant Decays 

Non-resonant D decays can enter into our signal modes with the same final particles. 
For example, non-resonant D + — > 7r + (7r + 7r~) can appear in the D + — > n + K s ,K s — > it + it~ 
mode. Also, non- resonant D + — > 7r + (7r + 7r - ^) can appear in the D + — > 7r + r)',r)' — > 7r + 7r _ ?7 
mode. To understand the backgrounds from non-resonant _D° or D + decays, we look at M^ c 
distributions in the invariant mass sideband regions of the intermediate resonances (K s or 
r)'). For decays, we follow the same procedure, replacing M^ c with M(D S ). The scaling 
factor, from sideband to signal region, is taken to be unity, as indicated by Monte Carlo 
studies. 

For the D° — * K S K S (or D° — > K s r)') mode, the scatter plot of K s candidate invariant 
mass against the other K s (or rj') candidate invariant mass is used to define a signal region 
and two kinds of sideband regions to remove the non-resonant decay background. Again, 
the scaling factor, from sideband to signal region, is taken to be unity. 

E. Systematic Uncertainties 

We have considered several sources of systematic uncertainty. Some are correlated among 
different decay modes. These include: 

1. the uncertainty associated with the efficiency for finding a track - 0.3% per track [2]; 

2. an additional 0.6% per kaon track is added [2], uncorrelated with item 1; 

3. the uncertainty in charged pion identification is 0.3% per 7r =l= [2]; 

4. the uncertainty in charged kaon identification is 0.3% per K ± [2], uncorrelated with 
item 3; 

5. the relative systematic uncertainties for 7r°, K s , and rj finding efficiencies are 2.0%, 
1.8% [2], and 4.0%, independent of one another, and independent of the first four- 
mentioned uncertainties; 

6. finally, among the correlated systematic uncertainties, there are the uncertainties in 
the input branching fractions of the normalization modes, 2.0% for D° — > K~ir + [2], 
2.2% for D + -> K-n + n + [2], and 5.8% for D+ -> K° S K + [3]. 

Note that for K s , with Kg — > 7r + 7r~, item 1 applies, as the tracks must be found, but item 
3 does not apply, as pion identification is not required for K$ — > 7r + 7r~ . 

The systematic uncertainties that are uncorrelated among the decay modes include those 
due to choice of signal shape and background shape. They range from ±0.05% for the 
cleaner decay modes to ±4.55% for the modes with substantial background. 

In the Table II we separately list, for each decay mode, the quadratic sum of the systematic 
errors excluding that from the normalization mode, and the error from the uncertainty in 
the normalization mode. 
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F. CP Asymmetries 

The Standard Model predicts that direct CP violation in D decays, e.g., a difference in 
the branching fractions for £>+ — > K + r] and D~ — > K~rj, will be vanishingly small. We have 
separate yields and efficiencies for D and D events, so it is possible to compute asymmetries 
Acp = (£>+ — B-)/(B + + B-), which are sensitive to direct CP violation in D decays. 
All systematic uncertainties cancel in this ratio, with the exception of charged pion and 
kaon tracking and particle identification efficiencies. Here the relative factor is the charge 
dependence of the efficiencies in data and Monte Carlo simulations [2] . 

For D° vs. D°, the only asymmetry we can measure is K~n + vs. K + tt~. That difference 
will contain a component from the difference in the doubly- Cab ibbo-suppressed decays D° — > 
K + n~ vs. D° — > K~ir + , as well as the component from the favored decays D° — > K~tt + vs. 
D° — > K + n~ . Our measurement does not separate these two possible asymmetries. 

VI. SUMMARY 

The obtained branching ratios, branching fractions, and CP asymmetries for all D — > PP 
modes are shown in Table II. The values we obtained are consistent with the world aver- 
ages [15] and for the suppressed modes, of better accuracy. No significant CP asymmetries 
are observed. 
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TABLE II: Ratios of branching fractions to the corresponding normalization modes D — ► K~ir + , 
D + — ► i^~7r + 7r + , and D+ — ► branching fractions results from this analysis; and charge 

asymmetries Acp- Uncertainties are statistical error, systematic error, and the error from the 
input branching fractions of normalization modes. (For D°, the normalization mode is the sum of 
D° -► K-TT+ and D° -► i^+vr" - see text.) 



Mode 




^mode /^Normalization (%) 


This result B (%) 


Acp (%) 


D° 




K+K- 


10.41 ± 0.11 ± 0.11 


0.407 ± 0.004 ± 0.004 ± 0.008 




D° 




K°sK°S 


0.41 ± 0.04 ± 0.02 


0.0160 ± 0.0017 ± 0.0008 ± 0.0003 




D° 




7T+7T 


3.70 ± 0.06 ± 0.09 


0.145 ± 0.002 ± 0.004 ± 0.003 




D° 




7r°7r° 


2.06 ± 0.07 ± 0.10 


0.081 ± 0.003 ± 0.004 ± 0.002 




D° 




K--K+ 


100 


3.9058 external input [2] 


0.5 ± 0.4 ± 0.9 


D° 




K° s 7r° 
K° s rj 


30.4 ± 0.3 ± 0.9 


1.19 ± 0.01 ± 0.04 ± 0.02 




D° 




12.3 ± 0.3 ± 0.7 


0.481 ± 0.011 ± 0.026 ± 0.010 




D° 




7r°rj 


1.74 ± 0.15 ± 0.11 


0.068 ± 0.006 ± 0.004 ± 0.001 




D° 




K° s i 


24.3 ± 0.8 ± 1.1 


0.95 ± 0.03 ± 0.04 ± 0.02 




D° 




n°r/ 


2.3 ± 0.3 ± 0.2 


0.091 ± 0.011 ± 0.006 ± 0.002 




D° 




7777 


4.3 ± 0.3 ± 0.4 


0.167 ± 0.011 ± 0.014 ± 0.003 




D° 




rjrj' 


2.7 ± 0.6 ± 0.3 


0.105 ± 0.024 ± 0.010 ± 0.002 




D+ 






100 


9.1400 external input [2] 


-0.1 ± 0.4 ± 0.9 


D+ 






3.35 ± 0.06 ± 0.07 


0.306 ± 0.005 ± 0.007 ± 0.007 


-0.2 ± 1.5 ± 0.9 


D+ 






1.29 ± 0.04 ± 0.05 


0.118 ± 0.003 ± 0.005 ± 0.003 


2.9 ± 2.9 ± 0.3 


D+ 






16.82 ± 0.12 ± 0.37 


1.537 ± 0.011 ± 0.034 ± 0.033 


-1.3 ± 0.7 ± 0.3 


D+ 






0.19 ± 0.02 ± 0.01 


0.0172 ± 0.0018 ± 0.0006 ± 0.0004 


-3.5 ± 10.7 ± 0.9 


D + 




K+T] 


< 0.14 (90% C.L.) 


< 0.013 (90% C.L.) 




D+ 




■K + r] 


3.87 ± 0.09 ± 0.19 


0.354 ± 0.008 ± 0.018 ± 0.008 


-2.0 ± 2.3 ± 0.3 


D+ 




K+T]' 


< 0.20 (90% C.L.) 


< 0.018 (90% C.L.) 




D+ 




ir + r]' 


5.12 ± 0.17 ± 0.25 


0.468 ± 0.016 ± 0.023 ± 0.010 


-4.0 ± 3.4 ± 0.3 


Df 




K° S K+ 


100 


1.4900 external input [3] 


4.7 ± 1.8 ± 0.9 


Df 




ir + ir° 


< 2.3 (90% C.L.) 


< 0.037 (90% C.L.) 




Df 




K° s tt + 


8.5 ± 0.7 ± 0.2 


0.126 ± 0.011 ± 0.003 ± 0.007 


16.3 ± 7.3 ± 0.3 


Df 




K+TT° 


4.2 ± 1.4 ± 0.2 


0.062 ± 0.022 ± 0.004 ± 0.004 


-26.6 ± 23.8 ± 0.9 


Df 




K+T] 


11.8 ± 2.2 ± 0.6 


0.176 ± 0.033 ± 0.009 ± 0.010 


9.3 ± 15.2 ± 0.9 


Df 




■K + r] 


123.6 ± 4.3 ± 6.2 


1.84 ± 0.06 ± 0.09 ± 0.11 


-4.6 ± 2.9 ± 0.3 


Df 






11.8 ± 3.6 ± 0.6 


0.18 ± 0.05 ± 0.01 ± 0.01 


6.0 ± 18.9 ± 0.9 


Df 




ir + r]' 


265.4 ± 8.8 ± 13.9 


3.95 ± 0.13 ± 0.21 ± 0.23 


-6.1 ± 3.0 ± 0.3 
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